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A COMPARISON OF COUPLING EFFICIENCIES FOR A STIX COIL AND AN m = 1 COIL 

by Donald R. Sigman 
Lewis Research Center 


SUMMARY. 

This theoretical and experimental study compares the ion-cyclotron wave generating 
characteristics of a Stix coil (which generates waves with azimuthal mode number m = 0) 
with those of a coil which produces primarily m = ±1 ion-cyclotron modes. The theo- 
retical work of J. E. Hipp, which predicted very good coupling for the m = 1 coil, has 
been extended to determine the scaling laws for plasma column radius Tp and coil 
wavelength A . Experimentally, an m = 1 coil '(X Q = 32 cm) and an m = 0 coil (X Q = 

40 cm) have been used to generate ion-cyclotron waves on a beam generated plasma 

12 _** 

column with electron density n g = 1x10 centimeter . Coupling resonances with peak 
efficiencies of approximately 40 to 50 percent were measured for both coils in low power 
(~10-kW) experiments. For equal power transfer to the plasma, the m = 0 coil voltage 
was more than a factor of two greater than that for the m = 1 coil. 

INTRODUCTION 

The use of a Stix (ref. 1) coil to couple energy to the ion-cyclotron wave (ICW) in a 
plasma has been studied both experimentally (ref. 2 and private communication from 

C. C. Swett) and theoretically (refs. 1 to 5) . This work has uncovered several features 

**. 

regarding the efficiency of coupling energy to the ICW and to the manner in which the. 
energy carried by the wave is distributed across the plasma column. It is well known 

(refs. 2 and 3) that the optimum coupling efficiency occurs at electron densities n in 

12 -3 e 

the neighborhood of 2 - 5x10 centimeter . For a fixed coil geometry and current, the 

energy coupled to the ICW varies approximately as the plasma radius to the fourth 

power. Because of the azimuthal symmetry of the Stix coil, the fields which heat the 

ions and the power flowing in the wave are largest near the outer edge of the plasma. 

Experiments (private communication from C. C. Swett) have verified that energy is 

coupled into the plasma by the Stix coil at an efficiency greater than 70 percent and that 

the wave generated in the plasma column has the dispersion characteristics predicted by 



the theory (ref. 3). 

Assuming the ICW to vary as e i(k z+m ® -u, t) j one can describe the wave generating 
current sheet of a Stix coil as a superposition of two traveling waves (one in the +z- 
direction and one in the -z-direction) each having m = 0. That is, 


. * 
J 



sin k ze 
o 


-icut 




\ 


( 1 ) 


(Symbols are defined in the appendix). Hipp, Kristiansen, and Hagler (ref. 6) have cal- 
culated the coupling of energy to the plasma for current sheets with m = 0, ±1. When 
m = -1 the waves travel azimuthally in the same direction that ions gyrate in the mag- 
netic field, likewise in the direction of electron gyration when m = +1. For m = ±1 the 
wave energy is maximum near the center of the plasma column and would thus be expec- 
ted to heat the innermost ions when damped in a "magnetic beach." Hipp presents the 
general equations for coupling to any m number mode, but only makes calculations fora 
particular situation (coil radius R = 3 cm; coil wavelength A = 277 /k = 15 cm; plasma 

° ° if _3 

column radius r = 2. 5 cm). He finds, that over the density range 1x10 centimeter 

14 P _3 

< n e <1x10 centimeter , the m = -1 current sheet couples to the ICW about an or- 
der of magnitude more power per unit current than does the m = 0 (i.e. , the plasma 
resistivity R L is larger for m = -1). Except at very low densities, the efficiency for 
coupling energy to m = +1 modes is much less than that for m = 0 or m = -1 modes. 

The coupling of energy to the ICW in a hot plasma by a short wavelength coil has 
been shown (ref. 5) to be inefficient because of the occurrence of ion-cyclotron damping 
underneath the coil. Therefore, in this report the calculations of Hipp (who assumed 
A =15 cm) are extended to calculate the resistivity R T for longer wavelengths. Also, 
the variation of R^ with plasma radius is determined. Experimentally, a Stix -like coil 
(i.e. , having a standing wave current sheet in the z-direction) with primary aximuthal 
numbers of m = ±1 has been constructed. Both this coil and a regular Stix coil are 
used to couple energy to the ICW in a beam plasma source. A comparison of the coup- 
ling efficiencies and radiofrequency coil currents and voltages is made. 


THEORY 


Geometry and Coupling Equations . 


The wave generating current sheet is assumed to be two wavelengths long and the 

i(k z+m0-wt) 

current density (i. e. , current per unit length) varies as e (see fig. 1). 
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/— Rf current sheet with radius, R; length, 2)^; 
/ and current varying as exp[i(k 0 z + m8 - wt)], 

/ V “ 2n/k 0 


L 



/ 

/ 

/ 


Plasma column with radius r^- 1 

Figure 1. - Radiofrequency (rf) coil and plasma model, cs-60950 ' 


The power coupled to the ICW in a cold collisionless plasma is given by Hipp, Kristian- 
sen, and Hagler (ref. 6) as 

+ 00 

0 \ ' r. sin 2 (k - k )X . f 

P av = 2lTiRl ) A 2 — Jdet[N ? ]K m (TR)cos xjs 

L-l W <k 0 -k n ) 2 Adet[A7 

n*0 


lmy 


d e t [N 7 ]K m ( T R)- 




det [ N 8 ] K m (TR ) 


sin i/A 


( 2 ) 


k=k 


The values k, for which det [A] = 0, are the wave numbers k n , of the natural modes; 
and the contour integral technique of Stix (ref. 1) was used to arrive at equation (2). The 
angle i p represents the direction of the current relative to the z -direction (t p - n/2 for 
Stix -like coils). All other terms are as defined by Hipp. The plasma resistivity is‘de- 
fined as 


R t 


av 


2ttI 2 L 2 


(3) 


Numerical Calculations 

Using equation (3) (with P from eq. (2)) computer calculations have been made to 

a. V 

determine R^ for all possible combinations of three wavelengths (A = 32, 40, and 
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48 cm), two plasma radii (r =2.5 and 5.0 cm), and three densities (n g = 5X10 11 , 

5x10^^, and 5xl0^ cm'^). For these calculations, the current sheet radius R = 7. 5 
centimeters, the frequency f — 6 . 5 megahertz, and the pitch angle \[/ = n/2. 

As in all ICW calculations, there is maximum coupling when the magnetic field is 
such that the wavelength of the natural mode X is close to that of the wavelength of the 
current sheet X . Table I gives R T with m = -1 for the various parameters discussed 

O i-i 

previously. The gas is deuterium. In each case, the magnetic field is that needed for 

maximum coupling. Table II gives maximum R^ with m = 0 and X Q = 40 centimeters. 

The data for table n is from reference 3 and the gas in that calculation was hydrogen. 

The quantity v in the last column of each table gives the variation of resistivity with 

plasma radius according to the empirical relation R^ cc where calculated values of 

R at r = 2. 5 and 5. 0 centimeters were used to determine v. 

L p 


TABLE I. - PARAMETRIC COMPARISON OF PLASMA RESISTIVITY 


R l (m = -1 CURRENT SHEET; DEUTERIUM IONS) 


Wave- 

Plasma radius, r p , cm 

V 







V 

2 

5 

5. 

0 


cm 

Plasma 

Nondimen- 

Plasma 

Nondimen- 



resistivity, 

sional 

resistivity, 

sional 



R L’ 

frequency , 

R L’ 

frequency , 



ohms 

n a 

ohms 




Electron density, 

n e , 5X10 11 cm -3 


32 

2.61X10" 2 

0.989 

1.72X10" 1 

0.990 

2.72 

40 

3.79 

.983 

2.32 

.985 

2.62 

48 

4.78 

.977 

3.04 

.980 

2.67 

64 

6.48 

.959 

4.24 

.969 

2.71 

80 

7.88 

.946 

5.30 

.950 

2.75 

Electron density , n g , 5xl0 42 cm"^ 

32 

4.89X10' 2 

0.924 

2.37X10' 1 

0.920 

2.28 

40 

6.95 

.895 

3.41 

.890 

2.30 

48 

9.04 

.864 

4.25 

.860 

2.24 

Electron density , n g , 5xl0 43 cm'^ 

32 

3.73X10' 2 

- 0.73 

4.90X10' 1 

0.64 

3.72 

40 

1.31X10' 3 

.69 

3.15 

.62 

7.90 

48 

1.26X10' 4 

.60 

9.40X10 -2 

.60 

9.51 


a n = 2rf/Q ci = w/B ci . 
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TABLE H. - PARAMETRIC COMPARISON OF PLASMA RESISTIVITY 


R L (m = 0 CURRENT SHEET; HYDROGEN IONS) 


Wave- 

Length, 

A , 
o’ 

cm 

Plasma radius, r^, cm 

V 

2.5 

5.0 

Plasma 

resistivity, 

R L’ 

ohms 

Nondimen- 
sional 
frequency , 

Plasma 

resistivity, 

R L’ 

ohms 

Nondimen- 

sional 

frequency, 

R a 

11 -3 

Electron density, n g , 5x10 cm 

40 

1.30X10" 4 

0.95 

4.13X10" 3 

0.97 

4.98 

12 -3 

Electron density, n , 5x10 cm 

40 

2.34X10' 4 

0.81 

5.45X10" 3 

0.81 

4.55 

13 -3 

Electron density , n , 5x10 cm 

40 

_5 

9.57x10 0 

0.43 

-3 

1.95x10 J 

0.43 

4.35 


a D = 2»f/B ci = w/S2 ci . 


Discussion of Tables I and II 

It is difficult to discuss the results in broad general terms. For example, the den- 
sity for optimum coupling depends on both coil wavelength and plasma radius; and, al- 
though not shown in table I, on coil radius as well. One general feature is that at high 
densities there is a rapid falloff in Rj^. The density at which this falloff begins is X Q 
and r p dependent, being lower for larger X Q and/or smaller r p . For example, at 
X Q = 40 centimeters and r = 2.5 centimeters, R. is reduced by a factor of 53 in going 

from n = 5x10 to 5x10 centimeter . At X =40 centimeters and r =5.0 centi- 

e o p 

meters, this factor is only 1.1. And, finally at X Q = 32 centimeters and Tp = 5.0 cen- 
timeters, R t is increased by a factor of 2.1. In Hipp's calculations (X =15 cm, 

r = 2. 5 cm, R = 3.0 cm), the optimum density is above 10 centimeter" and the rapid 
P 14 -3 

falloff occurs at n g > 10 centimeter Thus, to avoid low values of at high den- 
sity, one must use a short wavelength coil. 

Ip O 

Near conditions of good coupling (n = 5x10 cm and r = 2.5 or 5.0 cm) the fol- 

P 

lowing features are rioted: (1) Coupling improves with increasing X . However, before 
the optimum is reached, the coil length 2 X q becomes too long to be practical (i.e. , 

«200 cm). (2) Coupling improves exponentially with increasing plasma radius with the 
mean value of exponent v being 2.3 (i.e. , R T ~ rV This is compared to v = 4.6 for 

Li P 
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the m = 0 coil at n = 5x10 centimeter . It should be noted that these values of v 
e 

are only meaningful near the optimum density conditions. Thus, in designing an experi- 
ment, one must be careful not to extrapolate the data of tables I and II too far. 


EXPERIMENT 

Apparatus 

Figure 2 is a schematic of the experimental setup. The conf. g magnetic field at 
the’ midplane was variable up to 0.45 T and the mirror ratio was 2 to 1. The distance 


Rf power input 



CS- 60949 H lm 


Figure 2. - Experimental arrangement for comparison of rf coils. 


between mirrors was 3.0 meters and the field was constant over the plasma volume to 

within ±1 percent for an axial distance of ±1.2 meters from the midplane. A dc beam 

plasma source was used which produced a plasma of radius r ~ 4 centimeters and 

12 -3 P 

electron density n g « 1x10 centimeter" (as determined by a movable Langmuir 
prGbe). When the beam current was set at 100 amperes, the beam voltage was 95 volts. 
The typical background hydrogen gas pressure (no plasma) was 2.3x10 torr. 

The radiofrequency coil (either an m = 1 or m = 0) was excited by a 
14-kilowatt transmitter with an output frequency of 6.0 megahertz. Figure 3(a) 
shows the two-wavelength Stix coil which was used and the circuitry used to match 
the input impedance of the coil to the output impedance of the transmitter (50 ohms). 
The coil radius R = 7. 5 centimeters and the wavelength \ = 40 centimeters. 

Figure 3(b) shows the m = 1 coil which was used and also its associated matching 
newwork. Again the coil radius was 7.5 centimeters. Because it was easiest to 
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O' 
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L Infinite plasma column 

Stix coil 
Surface 



Current-sheet representation 



(a) Conventional stix coil. 



m = 1 coil 


Surface 
current, j* 



Current-sheet representation 


C ^Matching network 




1 / 

Transmitter 

_j 






(b) New m = 1 coil. 

Figure 3. - RF coils and matching networks used in experiment. 
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construct the m = 1 coil as a 2 ^ wavelength coil, it was necessary to shorten the 
wavelength to 32 centimeters to make it fit the existing test section. The current dis- 
tribution for this latter coil is a standing wave azimuthally. Therefore it represents a 
superposition of two traveling waves, one with m = -1 (which should be efficient in 
coupling energy to the ICW) and one with m = +1 (which should not be efficient). A 
grounded Faraday strip shield was used with both coils. 

The inductance of the m = 1 coil was much less than that of the Stix coil (0.364 and 
2.54 pH, respectively). However, it was found that the same variable matching capaci- 
tors and Cg could be used by simply moving one lead. 

• Figure 4 shows typical measured resonance curves for the Stix coil and the m = 1 
coil, respectively. The data were always taken under exact impedance matching condi- 
tions and always after a sufficient time to reach thermal equilibrium in the system (30 to 
60 sec). Resonances measured soon after the system had been open to the atmosphere 
were usually poor (see solid curve in fig. 4(a)). However, after several high power runs 



(a) Total rf transmitter output power plotted against Q using 
Stix coil. 



1.04 1.02 1.00 . 98 . 96 . 94 . 92 . 90 . 88 - 86 

Nondimens ional frequency, Q ■ w/n c j 
(b) Total rf transmitter output power plotted against Q using 


m-lcoil. 

Figure 4. - Resonance absorption curves. 
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were made, the data resembled the remaining curves in figure 4 and did not change from 
run to run or even from day to day. The value of fi for maximum coupling agrees well 
with theory. At n g = 1x10 centimeter , theory predicts the maximum (for the X Q = 
40 cm Stix coil) to be at fl = 0.962. The measured value was 0.958. For the \ Q = 

32 centimeters and m = 1 coil, the theoretical and experimental values are 0.974 and 
0.968, respectively. These slight discrepancies might be just experimental error or 
they might be due to the fact that coupling radiofrequency power to the plasma increases 
the density thus driving the resonance peaks to lower values of Q. 

The coupling efficiency CE for power transfer to the plasma is given by 


CE 


P. 


in 


P in + P loss 


(4) 


where P in is the power transferred to the plasma and P^ oss represents the circuit 
losses. Experimentally, P^ oss is determined by running with no plasma ; and P^ ^ 

(= P^ n + P^ oss ) is the net forward power when plasma is present. Table III is a summary 
of the resonant coupling efficiencies achieved with both coils at a variety of operating 
levels. 


Comparison of Stix Coil and m = 1 Coil Results 


Using the data of table III, we can make a comparison of the two coils by calculating 
the coil currents and voltages necessary to make P in = - P i oss = 6 kilowatts. 

TABLE III . - SUMMARY OF m = 1 COIL 
EXPERIMENTAL DATA 


Coil 

Peak coil 
voltage , 
V 

peak’ 

kV 

Rf circuit 
power loss, 

p 

loss’ 

kW 

Total of 
transmitter 
output 
power 

P total , 
kW 

Coupling 

efficiency, 

CE 

m = l 

5 

2.0 

3.7 

46.0 

m = 1 

7 

4.2 

7.6 

44.8 

m = 1 

8 

5.8 

9.9 

41.4 

m = 1 

9 

7.8 

13.2 

40.9 

Stix 

22 

6.3 

13.0 

51. 5 
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For the Stix coil, we get Vp eak 

= 9.5 kilovolts and I 


20.8 kilovolts and I = 154 amperes; and for the 
rms 


m = 1 coil. V 


peak 


rms 


= 490 amperes. It is interesting to calculate 


from the aforementioned the current per unit length in each coil, 


• * m rms 


(5) 


where N is the number of turns per section (Ngj.^ = 4, = 1). We get j ^stix 

= 30.8 amperes per centimeter and = 15.3 amperes per centimeter. Note that 

for the m = 1 coil the proper current to use in computing j* is 490 amperes/2 
= 245 amperes since the current arriving at each section must split into two parts. 

To determine whether the values of P^ n for m = 0 and m = 1 are related according 
to theory, we take a measured value of P. n (m = 1 coil) and apply various theoretical 
correction factors to calculate the expected P^ n (Stix coil). There are five correction 
factors; they are the following: 

(1) Correction factor for R^: Using the available theoretical data we can calculate 
R L (r p = 4 cm) for m = 0 and 1 and for current sheets varying as e i ( kz+me )_ At 

n = 5X10 11 centimeter -3 , R T (m = 0)/R T (m = -1) = 0.0153; and at n = 5xi0 12 centi- 

c a Jj J-j G 

meter , R L (m = 0)/R L (m = -1) = 0.0149. Thus, we will choose the correction factor at 

n = lxlO* 2 centimeter -3 as 0.015. 
e 

(2) Coil turns correction factor: The number of turns per section in the Stix coil is 

2 

4 and in the m = 1 coil is 1. Since P. ~ N , the correction factor is 16. 

in ’ 

(3) Fourier amplitude correction factor: The axial current distribution in neither 
coil is sinusoidal as was assumed in the theory. Thus, the relative Fourier amplitude 
(evaluated at k = k Q ) has been evaluated for each coil and the ratio of these Fourier am- 
plitudes is 1.75. 

(4) Standing to traveling wave correction factor: The m = 1 coil does not represent 

a traveling wave in the azimuthal direction, but a standing wave. Thus only one-half the 

o 

amplitude (m = -1) is effective in coupling to the ICW. Since P^ n ~ (j*) the correction 
factor here is 4. 

(5) Relative total current correction factor: Finally, we must correct for the fact 

that the current per turn was not the same in each coil. This correction factor is 
2 2 

41 (m = 0) T (m = 1) and will be evaluated later. The factor of 4 comes from the fact 
that, at each turn in the m = 1 coil, the current splits two ways. 

Applying all these correction factors we get 
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P in (m = 0) = PJm = 1) 


H^(m = 0) 
R L (m = 1) 


2 

f N ^ m - - - PH (Fourier amplitude factor)^ 

LN(m = 1)J 


2 

x (Standing-traveling wave factor)^ [ ~ ^ m ~ : 01 1 ( 6 ) 

Ll(m = 1)/2J 


From table HI we take the following m = 1 values: Vp ea ^ = 8 kilovolts, I rms = 
412 amperes, and P in = 4.1 kilowatts. Now, using equation (1) we calculate the Stix 
coil power P in (m = 0) for a current of 165 amperes: 


P. n (m = 0) = 4.1 kilowatts(0.015)(4) 2 (1.75) 2 


165 

.412/2 


7 . 8 kilowatts 


The measured value of P. (m = 0) was 6.7 kilowatts. The inclusion of m = +1 modes in 

in 

the calculation does not explain the 15 percent discrepancy. Thus we see that the theory 
is quite useful in predicting m = 0 coupling if m = 1 results are known, or vice versa. 


CONCLUSIONS 

A theoretical parametric study of ion-cyclotron wave (ICW) coupling using an m = -1 
traveling wave coil structure shows the following features: 

1. There is a density for optimum coupling and also a density above which coupling 
falls off rapidly. Both densities depend strongly on coil wavelength A q and on plasma 
radius r being higher for shorter wavelengths and/or larger plasma radii. 

r 

2. Near the optimum coupling density, the plasma loading resistivity R^ (which is 
proportional to coupling) varies approximately as with u = 2.6. For an m = 0 Stix 
coil v ~ 4. 5. 

3. Also, near optimum coupling density R^ improves with \ Q until a maximum is 

reached. However, the total length of the coil (L = 2A ) at this maximum is probably 

11-3 ° 

impractical (A q = 80 cm for n g = 5x10 cm , Tp = 5.0 cm). 

4. For a given coil current density j*, R^(m = -1) is higher than RjJm = 0) for 
almost all values of A q , r , and n g . However, the practical problems of designing a 
coil for high current density may be more severe for the m = -1 coil than for the m = 0 
Stix coil. 

Experimentally, a coil which generates both m = ±1 traveling waves and an m = 0 
Stix coil have been used separately to couple energy to the ICW. The coupling efficien- 
cies for the m = 1 coil varied inversely with the power level from 41 percent to 46 per- 
cent. The coupling efficiency for the m = 0 Stix coil was 51 percent. An analysis of the 

11 



different geometry of the two coils has explained the lower efficiency of the m = 1 coil. 
The larger theoretical value of R L (m = 1) has been experimentally verified; but, in 
order to achieve higher efficiencies, a different geometrical design is needed, with more 
turns per section and with an azimuthally traveling current sheet (m = -1). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 15, 1972, 

764-74. 
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APPENDIX - SYMBOLS 


C l’ C 2 
CE 

f 

I 

I 


beam 


I 


mis 


i* 


K (TR) 
m ' 

k 

k o 

L 

m 

N 

N 7’ N 8 

n e 

P in 

P 

loss 

P total 

P 

R 

R t 


coil 

^peak 


defined in ref. 6 

rf matching network capacitances 
coupling efficiency, eq. (4) 
rf driving frequency, hz 
defined in ref. 6 
plasma source beam current 
rms coil current 
rf current density 
magnitude of j* 

1L. 

modified Bessell function of the second kind (m in order) 

ICW wave number 

rf coil wave number 

one -half rf coil length 

azimuthal mode number 

number of turns /section, rf coil 

defined in ref. 6 

electron density 

power transferred to the ICW 

rf circuit power loss 

total rf transmitter output power 

background neutral pressure (no plasma) 

rf coil radius 

plasma resistivity, eq. (3) 

plasma radius 

defined in ref. 6 

time 

rf coil voltage (rms) 
peak rf coil voltage 



z axial position coordinate 

y defined in ref. 6 

9 azimuthal position coordinate 

A ICW wavelength 

A q rf coil wavelength, L 

p o permitivity of free space 

v defined by eq. (4) 

\p pitch angle of rf coil current relative to magnetic axis 

ft nondimensional frequency, y/ft^ 

ft £ ^ ion-cyclotron frequency, radians 

y wave frequency, radians 
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